Effects of perturbing forces on the orbital stability of planetary 

systems 



L.Kiseleva-Eggleton^ 

IGPP, Lawrence Livermore National Laboratory, L-413, 7000 East Ave, Livermore, CA 

94550 

IkisselevaOigpp . ucllnl . org 

and 



E.Bois 

Observatoire de Bordeaux, 2 rue d'Observatoire, B.P.89, F-33270, Floirac, France 



boisQobserv .u-bordeaux . f r 



Received ; accepted 



^Also Dept. of Physics, University of California, Davis 



- 2 - 



ABSTRACT 

We consider dynamical effects of additional perturbative forces due to the non- 
point mass nature of stars and planets: effects such as quadrupolar distortion and 
tidal friction in the systems of exo-planets. It is shown that these forces should 
not be neglected while modelling the dynamics of planetary systems, especially 
taking into account the undefined real masses of the planets due to unknown 
orbital inclinations and the unsatisfactory application of Keplerian fits to the 
radial velocity data in multiple planetary systems. 

Subject headings: celestial mechanics, stellar dynamics - planetary systems - 
stars:individual {v Andromedae) 

1. Introduction 

About 50 extrasolar planets have been discovered so far, among them: at least three 
confirmed planetary systems around the stars v And with three Jupiter mass planets 
(e.g. Laughlin & Adams 1999, Rivera & Lissauer 2000); HD 83443 with two Saturn type 
companions (Mayor et al. 2000); and GJ 876 orbited by two resonant planets (Marcy et al. 
2001). In addition there are two substellar-mass companions orbiting HD 168443 (Marcy et 
al. 2001, Udry et al. 2001), a system of Earth-mass planets around pulsar PSR 1257+12 
(e.g. Konacki et al. 1998), as well as yet-to-be-confirmed systems such as a second planet 
around 55 Cnc (Jayawardhana et al. 2000), planets orbiting Lalande 21185 (Walker 1996), 
and three Earth-type planets around pulsar PSR 1828-11. No doubt the number of detected 
planetary systems will continue to increase rapidly due to the improvement of observational 
detection techniques and to new space missions (e.g. Kepler, COROT, FAME, SIM etc.) 
scheduled for the next few years. 
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The majority of members of new planetary systems have eccentric orbits. This 
may increase the dynamical interaction between components and make systems chaotic 
and potentially hierarchically unstable, ending up with an ejection of one or more (in 
systems with more than 2 planets component) to infinity. Other main factors which 
add to uncertainty of the dynamical state of the systems are (a) the unknown value of 
orbital inclination i, which allows us to know only the lower limit of planetary masses rrio 
from the function nio — nip sin i, and leaves the real mass nip as almost a free parameter 
(for attempts of i estimation see Gatewood et al. 2001, Pourbaix 2001), and (b) the 
unknown relative inclination between planetary orbital planes. Also, usually the orbits 
in planetary systems are not Keplerian because of mutual interaction between planets. 
Therefore standard Keplerian fits to velocity observations are strictly speaking not vahd. 
This point was addressed by Laughhn & Chambers (2001) who suggested a new four-stage 
procedure for the determination of dynamical parameters of multiple planetary systems. 
This procedure includes multiple-Keplerian fits using a semi-analytic scheme followed by a 
final self-consistent pohsh with N-body interations. This method should be a substantial 
improvement to existing fitting techniques. However, in cases when close approaches of two 
or more components are possible, one should take into account the non-point-mass nature 
of the bodies and, as a consequence, the following perturbations to Newtonian gravity: (a) 
the quadrupolar distortion (QD) of the bodies due to their mutual gravity; (b) the further 
quadrupolar distortion due intrinsic spin of the components; (c) tidal friction (TF); (d) 
General Relativity. In this paper we present a few examples showing the influence of two of 
these perturbations - (a) and (c)- on dynamical stability of systems with two planets. 
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2. The model 

In order to estimate the dynamical effects of quadrupolar distortion of interacting 
bodies due to their mutual distortion (QD) and to tidal friction (TF), based here on the 
near-equlibrium approximation, on the planetary systems treated as N-body systems with 
extra forces due to QD and TF in addition to the Newtonian gravity, we applied the 
following formulation for the force of one body on the other, developed by Kiseleva, 
Eggleton & Mikkola (1998) - hereinafter KEM98: 
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Here are the stellar or planetary radii, and ak are the dimensionless dissipation rates for 
the two bodies. Q is a version of the apsidal motion constant, a dimensionless measure of 
the distortability of the body (star or planet). KEM found that for a polytrope of index n 
in the range < n < 4.95, Q can be approximated by the interpolation formula 

3 / n\ 2.215 
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We apphed this model to the v And planetary system. Recently Jiang & Ip (2001) 
confirmed once more that the innermost planet docs not affect very much the dynamics of 
the middle and outer planets, and so we ignored it in our simulations. For calculation we 
used the regularized CHAIN method (Mikkola & Aarseth 1993) with perturbations. The 
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actual numerical integration of the equations of motion are carried out by a Bulirsh-Stoer 
integrator with a timestep accuracy of 10"^^. 

In our simulations we always used initial orbital parameters for orbital periods and 
eccentricities from the Lick Data (Butler et al. 1999, Rivera & Lissauer 2000): Pc = 242 
days, Cc = 0.23 for the middle planet C and = 1269 days, = 0.36 for the outermost 
planet D. We always started simulations with all three components positioned in the 
same plane at apoastrons of their orbits and with the two orbits out of phase with each 
other by QOdeg (so we did not use the observed values of uj). Orbital parameters for the 
V And planetary system are not very precisely defined and may vary rather significantly 
(particularly the eccentricities) during the lifetime of the system (see below). Note, however, 
that differences between values of orbital parameters given by different authors may not be 
explained by real changes during a short observational time of a few years. Any real and 
significant changes require a time scale of hundreds years. 

3. Results 

In our first set of simulations we use the nominal values of the planetary masses from 
the Lick Data: for sini = 1 mc = 1.98Af,/, = 4.11Mj. The adopted mass of the star 
was M = 1.3 Mq. The relative inclination between the two orbital planes was taken to be 
(j) — 1 deg (models with higher relative inclinations were also tested but we will discuss 
those results in our next paper). For this set of parameters the system appear to be 
hierarchically stable over at least 10'' years, despite both orbital eccentricities fluctuating 
almost quasi-periodically within significant ranges: e^"^ — e^*" ~ 0.5, e^"^ — e^**" ~ 0.15. 
These results are in good qualitive agreement with other simulations (Laughlin & Adams 
1999, Rivera & Lissauer 2000). Neither TF or QD or their combination changes notably 
(except very minor details) the orbital evolution of this system. 
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However, because of the unknown values of the inchnation of the orbital planes of the 
planets ic and i^, which can differ very much from each other, even the mass hierarchy of the 
V And planetary system is questionable. If one accepts the values of ic = 173.7 deg ±3.8 deg, 
id — 28.7 deg ±16.8 deg proposed by Pourbaix (2001) with a big uncertainty, this will lead 
to a very different system with component C being no longer a planet {rric ~ 18Mj), and 
rriii ~ 8.5Mj (so the mass hierarchy is now changed). In agreement with earlier stability 
analyses (Stepinski et al. 2000, Rivera & Lissauer 2000) the new system is highly unstable 
and disintegrates ejecting the component D within the first 1000 years. The addition of QD 
increases significantly the lifetime of the system, although it does not change the chaotic 
and unstable character of its dynamics. 

In order to study more systematically the effects of QD and TF (in this work we study 
these effect separately) on dynamical stability of planetary systems unstable in point mass 
approximation, we exchanged the masses of planets C and D in the v And system, leaving 
their intial eccentricities and orbital periods unchanged. Without QD and TF this system 
is unstable ejecting the less massive (m^ = 1.98Mj) component D within ~ 3 x 10^ years 
(top panel of Fig. 1). QD and TF may significantly change the dynamical evolution of the 
system. The scale of these changes depends very much on the chosen values of coefficients 
Qi and ai from Eq 2. The value of Qi is defined by the polytropic index of the body 
(Eq 3). For the presumably radiative star of 1.3 we assume n ~ 3 with corresponding 

Fig. 1. — Eccentricity of the outer planet in the v And planetary system with masses of 
planets C and D exchanged. The system is unstable without QD and TF (top panel) 
and with relatively weak QD and TF (second and fourth panels), although QD increases 
the life time of the system. QD with coefficients = 0.08 and Qc — Qa — 0.2 and TF 
with a — 10~^ make the system stable over the considered time interval (10^ years), with 
quasi-periodic fluctuations of (third and bottom panels) 
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0.02, and for planets we took Q ^ 0.2 {n ~ 1.5). The dynamics of the system during 
the first ~ 1.5 x 10^ years is significantly less chaotic than without QD and the planet D is 
not ejected to infinity until t ~ 5 x 10^ years. But the situation changes really dramatically 
when we increase the value of Q* to 0.08 (so assuming a more convective interior of the 
star). The third panel of Fig. 1 show that in this case the system displays hierarchically 
stable dynamical behaviour over at least 10® years, with quasi-periodic fluctuations of both 
orbital eccentricities similar to ones of original v And system. 

The dynamical effect of tidal friction depends on its coefficients ctj. Kiseleva et al. 
(1998) found that for stars similar to the ones in the A Tau triple system (Fekel & Tomkin 
1982) the most hkely a ~ 10~^. However, for planets a can be significantly larger. We 
tested our model with a = 10"^ and a = 10~^ for all 3 bodies. The results shown on the two 
lower panels of Fig. 1 are totally different, a = 10~^ does not seem to improve the stability 
of the system. However, for a — 10~^ the bottom panel presents once more a hierarchically 
stable system with quasi-periodic behaviour of its orbital parameters such as e^. Dynamical 
evolution of models with strong QD (Q* = 0.08) and with strong TF [a = 10~^) over 10® 
years look in this case remarkably similar, despite different dynamical properties of these 
perturbations: TF is a dissipative force with respect to the total orbital energy and QD is 
conservative. However, such a similarity does not appear in other cases (see below), and we 
suspect that over longer time the evolutionary patterns with QD and with TF will diverge. 

We also studed models of v And with sin i = 0.33 for both external planets C and D, 
so the mass hierarchy of the nominal system is preserved and Mc = 6.53Mj, = 13.56Mj. 
The results are shown on Fig. 2. TF and especially QD significantly increase the hfetime 
(upto correspondingly ~ 5 x 10^ and ~ 1 x 10^ yrs) of this unstable system. 
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4. Conclusions and prospects 

Our strongest conclusion is that the effects of perturbative forces such as quadrupolar 
distortion and tidal friction should not be neglected when investigating numerically the 
dynamical properties of extra-Solar planetary systems, especially when their hierarchical 
stability is questionable and even weak additional forces may change the qualitive character 
of their dynamics. In our examples QD always had a stabilizing effect, but we can not claim 
that it always works this way and more studies are needed. 

We did not consider here the contribution of General Relativity and intrinsic rotation 
of the bodies which under some circumstances may be important. In any case, because of 
problems with the reliable determination of orbital parameters which cannot be consider 
as Keplerian in N-body systems, it would be very useful to apply a good dynamical 
chaos indicator (Lyapunov-type exponent, for example), which can distinguish a long-term 
dynamical instability from relatively short-term integrations, given a good-size sample of 
possible initial parameters. One such indicator was suggested recently by Cincotta & Simo 
(2000). Our first attempts to apply it to extra-Solar planetary systems are very encouraging 
and we are going to discuss the results in our next paper (Gozdziewski, Bois, Maciejewski 
& Kiseleva-Eggleton, in preparation). 

The authors are grateful to the John Templeton Fondation for the grant which 
supported the publication of this paper in ApJ. LK-E thanks the Bordeaux Observatory for 

Fig. 2. — Eccentricities of planets C (left) and D (right) in v And with sinic = sini^ = 0.33. 
Without TF and QD the system is destroyed within the first 20000 years (top panel). TF 
(middle panel) and cspecilly QD (bottom panel) significantly increase the lifetime of the 
system. The system with QD actually disintegrates after 10^ yrs. 
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hospitality. 
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